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SUMMARY 


The nonlinear characteristic differential equations applicable to 
a quasi-one-dimensional unsteady channel flow with friction and heat 
transfer are linearized and integrated in functional form for the par- 
ticular study of small perttirbations from ideal shock- tube flows. If 
the equivalence of lansteady- and steady-flow bomidary layers is assumed, 
the problem of determining the perturbations in the unsteady flow reduces 
to an evaluation of the drag of a flat plate in the equivalent steady 
flow. 


For air at initially uniform temperatiire , the theory eval-uated with 
an equivalent steady-flow turb\ilent-bo\mdary-layer skin-friction coeffi- 
cient predicts that shock attenuation increases with distance and that 
average values of static pressure, velocity, density, and Mach number 
at a fixed position in the hot gas increase with time, whereas average 
sonic speed simultaneously decreases with time at a fixed position. 

Experimental measurements of the shock attenuation with distance 
and static -pres sure variation with time at a fixed position for diaphragm 
pressiire ratios from approximately L to l8 gave good agreement with the 

theoretical predictions where a value of O.O 58 I X (Reynolds number) 
was used for the skin -friction coefficient. 


INTRODUCTION 


The shock tube has become a common aerodynamic testing facility 
because of its relative inexpensiveness and versatility. In a shock 
tube it is possible to obtain unsteady flows with a wide range of flow 
parameters, such as Reynolds number, Mach number, and temperature, that 
either co\ild not be obtained in steady-flow apparatus with the present 
temperature limit of known alloys or could be obtained only with massive 
and costly equipment. 

The various states of the flow of a perfect, nonviscid, nonconducting 
gas in a shock tube may easily be determined theoretically by application 
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of the basic equations of momentum continuity, and energy. (See, for 
example, refs. 1 to ) The theoretical states so determined are coni- 
cal in a distance-time sense (until the shock-tube end effects interfere) 
that is, parameters are invariant along rays which have the same ratios 
of distance from the diaphragm station to time elapsed since diaphragm 
burst. This perfect-gas flow is also characterized by two regions, of 
equal constant velocity and pressure but of different density and tem- 
perature, separated by a contact surface or entropy discontinuity. One 
of these regions has as its leading boundary the shock wave advancing 
down the tube into the low-pressure air at rest; whereas the other has 
as its trailing boundary the expansion wave progressing into the high- 
pressure section. Therefore, it is theoretically possible to obtain 
steady-flow aerodynamic data in the short -d\ir at ion steady flows that 
exist while either region is passing a fixed point along the shock tube. 

Unfortunately, the flow of a real gas in the shock tiibe departs 
significantly in many cases from the aforementioned theoretical flow. 

It is intuitively obvious that a real flow differs from the theoretical 
flow in that a pressure drop is required of a real fluid flowing in a 
tube and a further pressure change is required to account for heat trans- 
fer between the tube and the fluid. In addition, the flow in the region 
behind the entropy discontinuity has been found to be quite turbulent 
and erratic so as to make it of very limited value for most testing 
purposes. Although the flow between the shock wave and the entropy dis- 
continuity does not degenerate into large-scale turbulence, this flow 
is still affected by viscosity, heat transfer, other imperfect gas 
effects, and nonideal diaphragm burst. (See refs. 1, 5^ £^n^ ^•) Th^ 
nonideal condition of diaphragm burst may be minimized by the proper 
choice of diaphragm material for each particular initial set of condi- 
tions. However, the other effects are unavoidable, and it is to be 
expected that their magnitude will increase in importance in the high- 
temperature and high Mach number range, where the shock tube appears 
to be otherwise most advantageous. 

The attenuation in shock strength as the shock travels down the 
tube is the most obvious and easily measured deviation from perfect- 
fluid theory. This attenuation has been the subject of several experi- 
mental and theoretical studies (refs. 5> 5^ 6, and 7)* 

The theory of reference 5 ^7 t>e briefly outlined as follows. The 
unsteady flow of the hot gas between the shock and entropy discontin-uity 
is reduced to a quasi-steady flow by choosing a coordinate system fixed 
to the shock. The bo-undary-layer problem is then reduced to a laminar 
solution similar to the Blasius solution, except that the wall velocity 
is nonzero. The unsteady boundary condition of the receding entropy 
discontinuity is ignored as is the entire cold-flow region between the 
expansion and entropy discontinuity. The skin-friction and heat-transfer 
effects obtained from this solution are then averaged across the ass-umed 
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one -dimensional flow and the resulting changes in average momentum and 
average energy require the generation of moving waves. These waves are 
assumed to be generated at the entropy discontinuity and to travel 
unchanged in strength until they overtake the shock. The shock attenua- 
tion is then determined from the waves which have overtaken it since 
flow was first initiated. 

The theories of references 6 and 7 much less elaborate than 
that of reference 3 and are based on the assumption that the mass flow 
through the shock at a certain time is the same as the mass flow at the 
entropy discontinuity evaluated at the same time. The mass flow at the 

entropy discontinuity is detennined Trom the boimdajry -layer displacemen-t 
thickness and free-stream conditions corresponding to an unattenuated 
shock. This boundary-layer displacement thickness is detennined in ref- 
erences 6 and 7 fnom an extension, to a circular tube and rectangular tube, 
respectively, of the Rayleigh problem of the instantaneous acceleration 
to constant velocity of a flat plate in a fluid at rest. 

In an unpublished analysis made at the Langley Aeronautical Laboratory, 
Messrs. Paul W. Huber, Donald R. McFarland, and Philip Levine collaborated 
on a theory for shock attenuation based on the mass -flow decrement due to 
displacement thickness at the entropy discontinuity. In order to deter- 
mine the shock attenuation at a given time, the displacement thickness 
was evaluated in this theoiy for a time smaller than this given time by 
an interval equal to the wave-travel time from the entropy discontinuity 
to the shock. This displacement thickness was determined from empirical 
flat -plate steady-flow data at a distance from the leading edge of the 
flat plate equal to the distance from the diaphragm station to the entropy 
discontinuity. 

None of the aforementioned theories have given good agreement with 
experimentally detennined shock attenuation. Reference 8 reports the 
use of a chrono -interferometer to obtain the timewise variation of density 
in the two theoretically constant density regions. A density rise with 
time was noted for both regions. Comparison of the experimental data 
with the theory of reference 3 gave poor correlation. By employing vari- 
ous simplified flow models, the aforementioned theories have ignored the 
fact that unsteady waves are continually being generated by the effects 
of friction and heat transfer in the entire flow region. The result of 
these waves overtaking the shock is attenuation, and their motion along 
the tube results in variations in pressure, density, and velocity at 
given points. Therefore, in order to obtain a satisfactory understanding 
of the flow of a real gas in a shock tube, a theory must recognize and 
treat the wave system in the entire flow field. Such a theory has been 
derived in simplified form, and this theory, together with experimental 
correlation obtained at the Gas Dynamics Branch of the Langley Laboratory, 
is presented in this paper. 
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SYMBOLS 


A 

a 

B 


Cv 


D 

D( ) 

Dt 


area of flow 

velocity of sound 

constant defined by equation (dU) 

local skin-friction coefficient, 2r/pU^ 


coefficient of specific heat at constant volmne 
coefficient of specific heat at constant pressure 

if X Area 


hydrafilic diameter. 


Perimeter 


convective derivative, - + U ^ 

ot ox 


Hf 


heat added per unit mass by frictional dissipation 


h 
J 
I 

M = U/a 
Ma = U/^a 
Mp = U/ap 

Ms = Us/^ 

Mv = v/a<„ 

Mtf = w/ae 

N defined in equation (55) 

n reciprocal of velocity exponent in boundary layer, — = 


wall heat-transfer coefficient 

heat added per unit mass due to heat transfer and heat sources 
fixed distance along shock-tube axis 
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p 


p' 


p 




Q 


Q’ 


R 

Re 

r 

S 

T 

t 

U 

U* 


Us 

V 

w 

X 

y 


characteristic parameter, — a + U 

R 

effective characteristic wave parameter, p - — ^ 

7 R 

Prandtl nuinber 
static pressure 

static pressure immediately behind shock 

characteristic parameter^ a - U 

a ^ 

effective characteristic wave parameter, Q - 

7 R 

gas constant 
Reynolds number 
recovery factor 
entropy 
temperature 
time 

free-stream velocity 
shearing-stress velocity, ^r/p 
shock velocity 

velocity in boundary layer at y 

velocity of ray in x,t plot of figure 1, x^/t^ 

velocity of ray in x,t plot of figure 1, x^/t^ 

distance along shock tube from diaphragm station 
distance from surface 
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7 ratio of specific heats, Cp^Cy; assumed = 1.40 for con^jutations 

6 houndary-layer thickness; also Indicates differential quantity 

I characteristic derivative, ■■■ + (U ± a) y — - 

6t ot ox 



9 difference between wall temperature and adiabatic recovery 

temperatvure of flow 

|i coefficient of viscosity 

V coefficient of kinematic viscosity 

I distance flow has progressed along surface 

p density 



ap 


T wall shearing 


stress 

influence coefficients, defined in equa- 
tions (12), (13), and (26) 


Subscripts: 

Subscripts not included on the symbols defined above refer, in 
general, to values at points or within regions shown in figure 1. 
Exceptions to be noted, however, are as follows; 


o 

perfect-fluid value 

t 

at time t 

X 

at distance x 

exp 

experimental 

ref 

reference value 

theor 

theoretical 
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THEORY 


The basic steps followed in deriving the theory are first outlined. 
The characteristic differential equations for a quasi-one-dimensional 
unsteady channel flow are linearized on the basis of the perfect-fluid 
flow in a shock tube. The terms of these equations for a particular 
lamina of fluid in the shock tube are expressed as functions of the 
distance the lamina has moved along the shock tube. A transformation 
of variables of integration is then made whereby the time integral of 
the characteristic equation is replaced by a distance integral along 
the particle path. The values of the flow variables, such as pressure, 
velocity, temperature, and so forth, may then be found from the values 
of the integrated characteristic parameters. 

One method of evaluating the resulting integrals is to assume the 
equivalency of steady and unsteady flows based on particle-flow time. 
Application of this method reduces the solution for the characteristic 
equations to the simple confutation of the skin-friction integral (total 
drag) of a flat plate in steady flow. 


The characteristic equations for a quasi-one-dimensional unsteady 
channel flow as derived in appendix A are: 

^ = -a ^ ^ . a ^ R y - 1 ^ R ^ 

6t Dt 7 6t 7 Dt D 

^ ^ iQge ^ a ^ R y - 1 ^ R ^^^^f 

6t Dt 7 6t 7 Dt D 




Furthermore, under the assumption that the relation between skin 
friction and heat transfer is the same for both steady and unsteady 
flows, the convective derivative of entropy may be expressed (see 
appendix B) as 



Dt 



e p -2M27 

T r JD 


Ucf 


(5) 


For application to the study of flows in a shock tvihe , the preceding 
equations are linearized by assming that all coefficients and differential 
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operators on the right-hand side are constant at their respective 
perfect-fluid values. In other words, the lines making up the perfect- 
fluid characteristics network are used as paths of integration, and the 
various derivatives and coefficients are evaluated on the basis of the 
velocities, temperatures, and so forth, of the perfect fluid. Equa- 
tions (l) to ( 5 ) then become, if linearized on this basis (subscript o 
refers to perfect fluid) , 


6t 


/D loge a\ 




L 


(la) 


5t 


/d loge A 



+ 



2Uo^Cf 

D 


(2a) 




+ 


1 9o 

7 - 1 To 



Vf 


(5a) 


where 






1 a 


° 8x 


It is further assumed that, although the velocity of the leading 
edge of the expansion fan is correctly used as -a^, the fluid and 
characteristic velocities may be approximated by Uq, and (u ± a)^^ in 
the entire region from the leading edge of the expansion wave to the 
entropy discontinuity. This assumption, which greatly simplifies the 
computations, introduces errors that increase in magnitude as the 
expansion-fan region increases with shock pressure ratio. 

Equations (la) and (2a) are expressions for the derivatives along 
characteristics of Slope ^ = (U ± a)^ and equation (5a) evaluates a 
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derivative along a characteristic of slope Uq (particle path). If 
the length | is introduced and defined as the distance that a particle 
in the perfect fluid has traveled since its acceleration from rest by 
either the expansion wave or shock wave, the terms of equations (la) 
to (53’) Dia’y be evaluated as functions of | alone in either region a 
or (See fig. 1.) It should be noted that the functional relation- 

ship may be different for the different regions. When these assumptions, 
together with the definition of |, are applied to the region a between 
the expansion zone and the entropy discontinuity, the following relations 
are obtained from the plot of distance against time of figure 1: 


= -agt^ + .Uo(t - ti) 


i = Uo(t - ti) 


Substituting equation (1+) into equation (5) gives 


Uq 

i = — (x + agt) 

Uo + ^6 


Therefore 


^( ) _ d( ) _ Up d( ) 

hx. d^ ^ Uq + ag d| 


and 

M ^( ) 

^t d| ht Uq + ag d| 


(5) 


( 6 ) 


(?a) 


(Tb) 


Equations {"Jb.) and (7b) may be used to obtain the convective derivative 
as follows: 


at 



Un 


d( ) 


( 8 ) 
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and since the characteristic derivative is expressed as 

5t Dt dx 


( 9 ) 


this derivative may be written, for region a of figure 1, by substi- 
tuting equations (T^l) and (8) into equation (9) as 


5( ) 


= Ur 


iU + 


Ur 


6t d| 

6( ) _ Mg + T1 t 1 


Uq + 9-e 




< i( ) 


6t 


Mg + T 


Ur 


iLi 


(10) 


Similarly, for region p, which lies between the entropy dis- 
continuity and the shock wave , 


6( ) ^ CTMs - Mp + 1 d( ) 


6t 




M, 


P 


d| 


(11) 


When equations (3a), (8), (lO), and (ll) are substituted in equa- 
tions (la) and (2a) and the area term is neglected because of the 
averaging process described in appendix A, the following equations are 
found for regions a and p, respectively. In order to avoid needless 
repetition, the solutions for both 5P/5t and 6Q/5t are presented in 
one equation with both pltis and minus signs indicated for certain terms 
on the right-hand side. The upper signs apply to the characteristics 
with slope U + a and the lower signs to those with slope U - a. Thus, 


_1^ 

ae 




5t 

2 1 


D Ti ' 

.5t , 

V 

a 


jir - l)Mg T ^ Mg 


+ Mg2 + 

Mg + n 


0 a p - 2/5 

t" ^ 

■^a 


1 + 


Mg + T) i 1 


(Ma + Ti)(7 - l) 


Ucf 


a 


0p, 


a 




( 12 ) 
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6P 

f 

6t 



D ag 



S't. 

P 


- i)Mp ? i]Mp . * 


®P p -2/5 

0P 


aMg - Mp ? 1 


1 + 7 r 

(7 - 1) (aMs - Mpj 


Ucff 


D 


(15) 


These equations may be integrated along the characteristic lines. 
However, for ease of computation, since the slopes of these lines are 
assumed constant for each region due to the linearization process, inte- 
gration may proceed along particle paths by a sxii table change of variable. 
For example, from the geometry of the accompanying sketch of region a, 
the following equations apply: 


x(t) = -a^ti + U(t - tj^) 
x(t) = -agtjj + (U + ao,)(t - tb) 
x(t) = Uta + (U - ac,)(t - ta) 



Solutions of the first and second equations and of the first and 
third equations yield 


(t - ti) 


U + a^^ + a 
U + ag 


€ 


(t - tb) 
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and 


U + ap - a™ 

(t - ti) = (t - td) + 

U + a^ 

Therefore, from the definition of the length | in region a as 
I = U(t - ti), the following relations are obtained: 

Along ttb , 


I = (t - tb)U 


1 + M(X + Tj 


(Ih) 


Along 


+ 11 

U 1 + + Ti 


i = 


(t - td)U ^ + Utd 

Ma + n 


(15) 


(16) 


^ Mg + T] 

U Mq, + T] - 1 


(17) 


Therefore equation (l2) is integrated in region a with the change of 
variable specified in equations (l^l-) to (l7): 


Pt - Pb 




Ucf^bt 


or 

Pt - Pb ^Pg Ma+Ti ni(t) 

ae D + q + 1 Jo ^ 


(18) 
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and 


rfss) 5t 

a-e \^Va 


or 


Qt - Qd 


^Qa Hx + Tj , 

D Ma + n - 1 ^ 


(19) 


Similar relationships may be fo\xnd for region p. From the geometry of 
the accompanying sketch, it Is evident that x may be expressed as 
follows : 


x(t) = Ustj + U(t - tj) 
x(t) = Utf + (U + ap) (t - tf) 
x(t) = Ustk + (u - ap) (t - tk) 



Solving the above equations yields 


t 


ap(t - tf) 

“ U - U, + * 


and 



Ug - U + ap 
Us - U 


(t - tk) 


Consequently, since in region p, | = u ^t - tjj , 
apply: 


the following relations 
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Along tf t , 


I = 


gMg - Mp - 1 
aMg - Mp 


U^t - tf) + Utf 


( 20 ) 


d| - Mp 

O U = — 


U 0Mg - Mp - 1 


(21) 


and along y 


^ ^ 1 ^ aMg - Mp _ X 

0Mg - Mp 


( 22 ) 


d| - Mp 

U 1 + aMs - Mp 


(25) 


Therefore, equations (15) (20) to (25) can be used to obtain 


Pt - Pf _ 


I i$). “ ■ T r, 


,8t 


0Pg aMs - Mp nl(t) 


D aMg - Mp 


- 1 Jutf 


Cfp d| 


(24) 


and 


Qt " 


aMg - Mp 


6t = ^ r Ucf 

Jtv ^ Jtf 


,5t 


D aMg - Mp + 1 Jo 


i|(t) 


Cfp d| 


(25) 
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The entropy at aiQr point (x,t) in region a or p may be found sin 5 )ly 
by integration along a particle path (since ^ is the product of U 
multiplied by flow time, d| = U dt), so that from equation (3a) ^ 


^t ~ Sfa 

R 

St - Sj 
R 




27 

TT 


+ 




Cf d| 

P 



Cf d| 
a 

P 


(26) 


It is now necessary to consider the transmission and reflection 
of the characteristics at the entropy discontinuity separating the 
fluid which was initially at a high pressure behind the diaphragm from 
that which was initially at a low pressure ahead of the diaphragm. These 
effects are ln5)ortant, even in a linearized solution, because of the 
large differences in v which may exist in the fluids on each side of 
the discontinTiity. 


If the subscript o refers to conditions in a perfect fluid (that 
is, immediately after diaphragm burst in real fluid) and the notation 
of figure 1 is followed, it is shown in appendix C that, if A « 1 
and 7 is the same in both a and p. 


R 



R 


(1 + A) = 



-) (1 + X) 

/ o 


(27) 


Then, 

6P ^ ^f - Pc ^ ^ ^ % - Qc 
S^o PpQ - Puq 5Qo QPq - Qoo 


= 1 + 


- Pgp + Qf - Qp^ 


+ 





(28) 
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The second term on the right-hand side indicates the change in the 
characteristic value which would occur as a result of reflections and 
refractions even without any change in entropy. The third tenn evaluates 
the additional variation due to the changing entropy . 

The reflections at the shock wave of the incident waves are of 
higher order and may be ignored in this solution. 

The following equations result from the definition of P, 
and s/R and subsequent linearization with the limitation to small 
perturbations : 


~ ^ P + Q 

ae C ag 


(29) 


IL= ^ ~ ^ 

2ag 


(50) 


V_ 

Po 





a - ap 
7-1 ap 


S - Sp 
R 


+ 


« 1 


7 ^ ~ ^o ^ ~ ^ 

2 Spy ag 


S 


" ^o 


R 




(51) 
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_p_ 

Po 



S-So 

R 


* 1 


+ 


1 

2 ao\ 


Pq + Q - Qq\ 
^6 / 



(32) 


Consequently, the description of the flow at any desired distance 
time may be determined by combining equations (29) to (32) with 
equations (l8), (19), (24), (25), (26), and (28) with the proper limits 
of integration determined by geometry on the x,t plot. 


It is to be noted that the values of P, Q, and S/R must also 
be determined at the intersection of the entropy discontinuity with the 
characteristic passing through the point (v or w) at which the flow 
properties are to be determined. This step is required to compute the 
transmission and reflections of the characteristic parameters as deter- 
mined by equation (28) . 


In region a at the point Xy = Wty = (see fig. l), the 

equations become: 

p M^+1 

% + q Mq+t] 


Pw ” Pfa 


- n 


Pa D + T) + 1 Jo 


TjUty 


Cfa 


(33) 


Pa - P, 


d - ../■a = I 

ae a D M + Ti + 1 Uq 

a 


fl+TlMv/-Ma)utw 


Cf^dl 


(34) 


Qe - Qm ^e " 1 ^^s " 

Pi 


i+,,v«a)uv 


ag P aMg - Mp + 1 ^0 


Cfp d| (35) 




Qw - Qd 


d 1 ^ + n 


P + T1 - 1 
a * 


M^+l 




Ma+4 

(l+qM„-Mc)ut^ 




(36) 
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Sy - Sj 
R 







(57) 


Qd - Qe ^ 1 ^ (P(i - Pq) + (^ - %) ^ 

^ - Qpo *^Po 



^3o/^ 





(l+T)M^-M^)uiv 



d| 




\l+riVMa)uV 


a 


^0 


‘a 


<i| 


(58) 


It should be noted that the neglect of the expansion fan in the 
derivation is reflected in the term 1 + t)Mw - appearing in the 

limits of integration of the above equations. Values of such that 

the point w is in the expansion fan make the factor 1 + 
negative and give only meaningless answers. 


Equations similar to those for region a are obtained for region 3 
at Xy = Vty = My-aooty as follows; 


8c - Pb ^ M 1 ^ ^ 


l+M3-aMy)utv- 


Cf^dl 


Pv - Pf _ 


. 1 ctMs - M g_ 

D aMg - Mp - 1 


^ <TUty 

aMg-Mp ^ 

(l+Mp-aM^)uV ^ 


Qf - Qfi _ ^ 1 dMg - Mp P(l+Mp-aM^)uty 

a^ ^Qp D cjMg - Mp + 1 Jq 


(59) 


m 


Cfp d| (41) 
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Qy - Qk _ 


crMs - Mp 

•P 5 aMg Mp + 'l Jq 


aMg-Mp 


<TUty 


Cfp d| 


(42) 


Pf - Pc 
^Po " % 


= 1 


, (Pc - Fb) , 




(r - 1) ®Pc/^oo 




D 


/Yl+M3-aMy)ut^ 


cfp 


- 1 



=fa<15 


Sy - Sj Sy - S 


R 


= 1 






aMg-Mp 


<TUty 


Cfp d| 


(43) 


(44) 


The solution to the small-perturhation problem of shock-tube flow 
has now been reduced to the evaliiation of integrals of the form 

J* Cf (^)d| where | is simply a length formed by the product of the 

flow time multiplied by the fluid velocity. Any analytic or graphical 
method giving Cf as a function of | may be used in the evaluation 
of these integrals. 

In order to evaluate equations (55) to (44) , it is necessary to 
know the state of the boundary layer at a point in the shock tube which 
is denoted by a corresponding point |(x,t) of the characteristics 
diagram. The boundary -layer problem is a "hybrid" of the Rayleigh 
problem of the instantaneous acceleration of a flat plate and the Blaslus 
problem of the steady flow over a semi-infinite flat plate. It is 
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similar to the Rayleigh problem in that, at a given instant of time, all 
the flow between the shock and expansion wave has been accelerated from 
zero velocity to a constant velocity but is different from the Rayleigh 
problem in that the fluid has been in motion for varying time durations. 

It is similar to the Blasius problem in that, along a particle path in 
x,t coordinates, the boundary layer increases rearward from the shock 
or "leading edge," but is different in that the leading edge is not sta- 
tionary. A solution to this hybrid problem for laminar boundary layers 
was foTond in reference 5 for region p by reducing the unsteady-flow 
problem to a steady-flow problem by choosing a coordinate system fixed 
to the shock wave. The problem then became the Blasius problem modified 
for a nonzero wall velocity. In a similar manner a solution could be 
found for region a if the expansion wave were assumed to have zero 
thickness and the coordinate system were fixed to the leading edge of the 
expansion wave. 

For laminar flows, the solutions to the Rayleigh problem, the Blasius 
problem, and the modified Blasius problem of reference 3 all show the skin 
friction to be inversely proportional to the square root of the time that a 
fluid particle in the free stream outside the boimdary layer has been in 
motion over the plate. It should be noted that the constants of propor- 
tionality are different in each of the above cases. No solutions to the 
turbulent-boundary-layer equivalent of the Rayleigh or modified Blasius 
problem are known to the authors. One might assume, however, that the 
state of the boundary layer at a point in a turbulent flow (as well as 
in a laminar flow) over a flat plate with zero pressure gradient may be 
expressed as a function of the time interval that the fluid has been in 
motion over the flat plate or of the distance the outer fluid has moved 
along the plate. 

In the absence of a modified Blasius solution in region a and of 
any solution whatsoever for the analogous problem of tinrbulent-boimdary- 
layer flow in either a or p, an approximate method of evaluation for 
the integrals may be obtained by extending the above-mentioned time- 
dependency assumption further and assuming that the linearized unsteady 
flow in a shock tube has the same properties as an equivalent steady 
flow defined in the following manner: The properties of a lamina of 

fluid in the shock tube which has been in motion with a velocity Uq 
for a time t are equivalent to those of a lamina of fluid which has 
progressed rearward for a period of time t from the leading edge of 
a semi— inf ini te flat plate in a steady flow with free-stream velocity Uq. 
Under these conditions | becomes the distance from the leading edge of 
the flat plate. 

The introduction of this assumption reduces the small-perturbation 
solution to the singjle evaluation of the integral of the skin-friction 
distribution along a semi-infinite flat plate in steady flow. This 
integral may be evaluated if the relationship between cf and | is 
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prescribed either in analytic form or as a curve of cp plotted against 
I or of Cf against Re. The fonner method of relationship will 
usually permit evaluation of the integral in closed form, whereas the 
latter method may require graphical Integration. 

A theory has now been obtained for small -pertvirbation flow in a 
shock tube with the following assumptions: 

(1) Variations from perfect fluid flow shall be small . 

( 2 ) Averaged values are used for quasi -one -dimensional flow. 

( 3 ) Averaged values of p and U which identically satisfy the 
continuity equation are also assumed to satisfy the momentum equation. 

(4) The frictional dissipation is approximated as the product of 
the average velocity times the wall-shearing force. 

(5) The perfect-fluid characteristic net, velocities, densities, 
and so forth axe used to obtain linearized values. 

(6) The expansion fan is treated as a "negative shock wave" in that 
the final values of velocity, pressure, and so on, are assumed to exist 
immediately after the leading edge of the expansion fan. 

(7) If numerical evaluation of the theory is performed on the basis 
of equivalent-flat-plate steady flow, the following assumptions are also 
introduced: 

(a) The steady and unsteady skin frictions are equal for a 
fluid which has traveled a given distance over a flat plate in 
steady flow and for a fluid which has traveled the same distance 
along the shock tube in unsteady flow. 

(b) The boundary layer is small relative to the width or 
height of the shock tube. 


EXPERIMENTAL APPARATUS AND PROCEDURE 


Experiments to determine the validity of the theory were performed 

in a high-pressure shock tube 2 inches high by li inches wide in the 

Langley gas dynamics laboratory. Air at room temperature was used for 
both high- and low-pressure sections of the shock tube. Measurements 
of shock velocity and pressure-time variation were taken at various 
points along the shock tube. 
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The low-pressure side of the shock tube was at room conditions ^ 
whereas the high-pressure chamber was filled to pressures of 45^ 10, 

95, or 250 pounds per square inch gage. Either one or two sheets of 
soft brass shim stock O.OOI25 or 0.0015 inch thick was employed as a 
diaphragm. The number and thickness of the brass sheets were determined 
so that the diaphragm was near its breaking stress at the particular 
pressures of the test. Consequently, optimum diaphragm burst would 
result when the point of the rupturing mechanism pierced the center of 
the diaphragm. It has been found that Imperfect diaphragm bursts will 
result in both increased shock attenuation and scatter of experimental 
data. 


The shock velocity was measured at distances of 3»25^ 5*25^ 7*25^ 
9.18, and 11.18 feet from the diaphragm by inserting shock-tube sections 
of various length between the diaphragm and the glass -walled test sec- 
tion. A block diagram of the experimental apparatus is shown in fig- 
ure 2. Velocity was determined from measurements of the time required 
for the shock wave to traverse the 1.206 feet between the two li^t 
beams of the optical systems. (See figs. 2 and 5*) Each optical system 
consisted of a direct-current automobile lamp with line filament alined 
vertically at the focal point of a two— inch— diameter, seven— inch focal- 
length lens. The light emitted from the lens was masked to a vertical 
slit approximately O.050 inch wide and 1 inch hi^. Stray li^t on the 
upstream side of this beam was cut off with a sharp knife edge. The 
beam then passed perpendicularly through the test-section window and 
impinged on the edge of a second knife edge about ten inches from the 
test-section window. When the plane shock wave reached the upstream 
edge of each beam, li^t was refracted off the second knife edge onto a 
phototube. The resulting signal of each phototube was amplified by the 
circuit shown in figure 4. The output pulse from each thyratron (2D21) 
was used to trigger one channel of an eight-megacycle counter ^^chronograph. 
The pulse from the upstream system triggered the chronograph "start" 
circuit and that from the downstream system, the "stop" circuit. (See 
fig. 2. ) The time interval could thus be measiored to within 
±1/8 microsecond. 

Pres s\ire -time records were obtained for the four nominal pressure 
ratios at a distance of 8.15 feet from the diaphragm by employing a 
flush-mounted capacitor-type pressure pickup (Rutishauser Electronic 
Pickup Indicator, Type ST-127A) and associated apparatus, as shown 
schematically in figure 2. The signal from fluctuations on the pressure 
pickup was fed into the Y-axis of one beam of a four-channel cathode-ray 
oscillograph, and the resulting deflection of the beam was photographed 
by an NACA synchronous drum camera rotating at 1,800 rpm. A 1,000-cps 
gj^n0_\/ave timing trace was fed into another channel of the oscillograph 
in order to obtain a time base. 
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The oscillograph beams could be turned on and off electronically by 
application of a suitable signal to the Z-axis of the oscillograph. This 
signal was obtained from a piezoelectric crystal pickup, a preamplifier, 
and an NACA electronic shutter circuit. (See fig. 2. ) The crystal pickup 
was located 6.4 feet from the diaphragm so that the oscillograph would 
record the ambient pressure ("zero trace") before arrival of the shock 
wave. The time diiration of the oscillograph trace was adj-ustable up to 
20 milliseconds. 

A minimum of five test runs was made for each combination of nomi- 
nal pressure ratio and distance to the median position of the velocity- 
measuring system. The chronograph time interval, the shock-tube-wall 
temperature, and the barometric pressure were recorded for each test. 

Pressure-time records at x = 8.15 feet were taken, along with 
simultaneous velocity measurements at a median position of 9.25 feet. 

The aforementioned data were recorded and, in addition, a static cali- 
bration of the pressure gage was performed for each run. This calibra- 
tion was made by applying pressure to the shock tube in predetermined 
increments, and at each pressure level the oscillograph Z-axis was ener- 
gized so that the static-pressure cadibration trace of the gage could be 
photographed to place a press\ire scale on the drum-camera-film record. 


The experimental shock strengths shown in figure 5 were determined 
from the time interval At seconds for the shock to pass the slits of 
the measuring station a distance AZ feet apart as follows: 


RESULTS AND DISCUSSION 


Evaluation of Experimental Data 




which becomes, for 7 = 1.40 and dry air. 



The pressure was plotted at a value of x halfway between the phototubes. 
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The experimental data required some sli^t adjustment before it could 
he presented in figure 5, each part of which represents a common diaphragm 
pressure ratio or theoretical perfect-fluid shock pressxire ratio. Since 
the tests pertinent to a given nominal diaphragm pressure ratio were made 
with the high-pressure chaniber at the same gage pressure, variations in 
atmospheric pressure from time to time caused a variation in diaphragm 
pressure ratio. The required adjustment was made in accordance with the 
relations that follow (notation of fig* l)« 

Since the equation relating shock pressure ratio to diaphragm pres- 
sure ratio (see refs. 5, k, and ?) may be expressed in functional form as 


Zi = f(^\ (kl) 

Poo \PqoJ^ 

then 



and 



adjusted to 
ref. pg/p„ 



(49) 


The reference value of Pg/Poo each nominal diaphragm pressxore 

ratio was taken as the value existing when the pressure-time records were 
obtained by the capacitor pickup. 


Since each of the velocity symbols represents an average (on a given 
day) of from five to ten test runs, in which the scatter in any set taken 
on the same day is less than the scatter among averages of data taken on 
different days, the reason for the discrepancy between averages at the 
same distance and same theoretical pressure ratio is unknora. Te^erature, 
viscosity, and humidity corrections, determined from the linear theopr 
and applied to these experimental points, are insufficient to the 

averages coincide. The cause of the apparently high experimental pre - 
sure ratio at x = 5-2 feet (fig. 5) is also unexplainable at present. 
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The pres sure -time traces of the capacitor pickup located 8.13 f^et 
from the diaphragm station are shown in figure 6. The linear- theory 
values computed at x = 8.0 feet are plotted (circled points) for com- 
parison at millisecond intervals after the arrival of the shock. The 
shock pressure ratio obtained from concurrent velocity data for each 
particular run falls on the initial theoretical point within the accuracy 
of plotting. The velocity data^ obtained at x = 9-23 feet^ were adjusted 
to 8.13 feet by drawing a curve through the velocity point at 9*23 feet 
parallel to the experimental curves of figures ^(a,) to 5(^)i thereby 
compensation was provided for any scatter in the particular test. 


Evaluation of Linear Theory 

Constants eng^loyed . - The equations of the linear theory are evaluated 
on the assumption that standard atmospheric conditions exist in the low- 
pressure chamber unless specified otherwise. Temperature recovery fac- 
tors r of 0.90 and O.85 are used for the turbulent and laminar boundary 
layers, respectively. A Prandtl nunber Pp of O.7I ^ value of the 

Blasius constant B of 8.7O are assumed. 

The temperature-viscosity relationship 

iL = (50) 

T + 225.2 \Too/ 


was expanded for computational convenience to give the approximation 



(51) 


which is considered valid in the temperature range near 520° Rankine. 


Relation between Cf(|) and VarioTis theoretical, empirical, 

and experimental relations are available for relating cf to ^ for 


evaluation of the integral 


/'cf(|)d|. 


If the assumption of an 


equivalent-flat-plate steady flow is en^jloyed and the boundary-layer 
profile obeys the law 



( 52 ) 
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it is shown in appendix D for an incompressible turbulent boundary layer 
and for n and B Independent of | that 


2 2 


n 

n+3 / ^ \n+3 

(n + 2)(n + 3)B^ 





Then, the integral may be evaluated in closed form as 



Cf d| 


' "(c) 


2 

n+3 


n + 3 
n + 1 


n+1 

,n+3 


^2 


( 5 ^) 


Similarly, if the equivalent steacSy flow is used with an incompressible 
laminar boundary layer (see ref. 9)f the solution is 


Cf(|) = 0.664 



1/2 


(55) 



(56) 


If the moving wall or modified Blasius solution of reference 3 is 
used, the expression for Cf(|) in region p reduces to 


Cf(|) = - 


Ai 


(57) 
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In whicli A2/A]_ are constants evaluated in reference 5* Then 


c,(i)di = 


2/^Ys/vf/2(7l/2l^^ 

\Ai/irwy li J ii 


(58) 


Since reference 3 does not consider region no equivalent expression 
for this region is available. However, if the same general assumptions 
were applied to region a as were applied to region 3, evaluations 
of Cf(|) would again yield an equation in which Cf(^) is proportional 



In order to determine a possible correct relationship between Cf(|) 

and I to be used in evaluating the linear theory for the range of 
experiments covered in this paper, the curves of theoretical and experi- 
mental shock pressure ratio against distance of figure 5(c) were used. 

In addition to the experimental points obtained for a nominal diaphragm 
pressure ratio of 7-^55 (shock pressure ratio 2.525)^ theoretical curves 
are plotted by assuming an equivalent steady flow for (a) laminar skin 

friction (cf = (b) turbilent skin friction (cf = 0.058lHe 

n = 7) with and without heat-transfer effects; and (c) an experimental 


curve (from fig. 88 of ref. lO) giving the integrated value of Cf, 

which includes the effect of transition of the boundary layer from laminar 
to turbiilent. Since, as mentioned previously, the modified laminar 
Blasius solution of reference 3 not evaluated in region a, no theo- 
retical curve can be drawn for the corresponding Cf(|). However, due 
to the laminar assumptions of reference 3^ "the attenuation in shock pres- 
sure ratio would be proportional to \fx if Cf(|) were evaluated in 


region a with similar assumptions to those used for region p. Con- 
sequently, in order to determine whether such a modified Blasius treat- 
ment could possibly yield the correct relation between cf and |, a 


curve of the type 




(Constant) Vx was fitted to the experi- 


mental data at x = 9 feet and the constant was found to be O.0626. 

This resulting empirical curve, — ^ = 2.5^3 - 0.0626\/x, is also plotted 

1^00 


in figure 5(c)- 
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It is evident from an inspection of figure 5(c) that best agreement 
between theory and experiment is obtained either from the txirbulent skin- 
friction law including heat-transfer effects or from the experimentally 
determined skin-friction curve of reference 10. The reason for the simi- 
larity of the shock pressmre curves obtained by these two methods is 

simply that, for < 20 x 10^, the turbulent skin-friction law 
^Cf = 0.058lRg“^/5^ for a value of n = 7 when integrated closely 
approximates the curve in reference 10 except at low values of Rg , 


where the laminar and transition effects are apparent. These latter 
effects resiilt in the inflection in figure 5(c) of the c\irve based on 
the drag curve from reference 10. Inasmuch as these two skin-friction 
curves resiilt in almost the same theoretical shock-pressvure attenuation, 
either could be used for further evaluation of the theory. However, 
since an analytic closed-form integral is more convenient for computation 


in this case, the skin-friction law Cf = O.O58I 



- 1/5 


was used for 


evaluation of the linear theory for the other figures in this paper. 


It is also evident from figure 5(c) that the equivalent steady- 
flow laminar -boundary-layer skin friction predicts shock attenuation 
much less than that meas\ired. Furthermore, the empirical curve (atten- 
uation proportional to - which might be said to represent the shape 

of laminar solutions in general, including the modified Blasius solution 
does not indicate the trend of measured shock attenuation. Thiis, it 
might be inferred from these conditions that the boundary layer is turbu- 
lent in most of the shock-tube -flow regions a and p at these values 
of Rg. 


Comparison of Theory and Experiment 

Shock attenuation with distance from diaphragm .- The theoretical 
shock pressure loss and shock attenuation coefficients are presented as 
functions of perfect-fluid shock pressure ratios in figure 7- These 
relations are used to predict the attenuation of the shock with distance 
traversed from the diaphragm station in figure 5 and the theoretical 
resiilts are found to be in good agreement with the experimental data. 
Altho-ugh the theoretical points are computed for standard atmospheric 
conditions, the experimental points are adjusted to the nonstandard atmos 
pheric conditions existing when the pressure-time runs were made. In 
order to determine the magnitude of error so introduced, the theory was 
evaluated for the nonstandard (t = 5^2°, p = l4.7 Ib/sq in.) conditions 
for the diaphragm pressure ratio of I7.915 where the error introduced 
would be largest. The adjusted theory indicated by the dashed ciirve of 
figure 5(d) shows a sli^t increase in attenuation which agrees somewhat 


NACA TN 5575 


29 


better with experiment than the "standard atmosphere" theory. The 
adjustment for nonstandard conditions would be even less for the other 
pressure ratios. (See fig. 7(t>).) 

Use of the linear theory with a value of Cf = 0. appears 

in figure 5 to overestimate the attenuation slightly at the lower shock 
pressure ratios (below 2.5) ^th the opposite effect apparent at hi^er 
ratios. This trend could be due to any of a number of simplifications 
made in the theory and computations. Two of the more obvious factors 
are the Reynolds number and compressibility effects. As the diaphragm 
pressure ratio increases , the Reynolds number of the flow per unit length 
increases. Since, for the computations of this analysis, Cf is assumed 

proportional to R^“ 2/(n+5) 

^hich is equal to R0”^/5 for n = 7^ 
and since steady-flow experiments in general indicate a trend of Cf 
proportional to smaller negative powers (larger n) with Re increasing^ 
there is the possibility that values of n should be Increasing some- 
what with Re and diaphragm pressure ratio. Compressibility would tend 
to reduce the ratio of the compressible cf to the incompressible Cf 

and consequently would tend to reduce the attenuation as the flow Mach 
numbers increased with diaphragm pressure ratio. 

The neglect of the spreading of the expansion fan must also be con- 
sidered, since this omission must introduce an error which increases 
with shock pressure ratio (expansion-fan size). 

Other possible factors open to question are the use of constant P^, 
the evaluation of heat transfer on the basis of v of the fluid in the 
center of the flow (for example, not at some intermediate temperature 
between the wall and stream), and the use of r = 0.90(0.85). Howevel*, 
because of the ass-umptions in the one -dimensional averaging process, in 
the linearization procedure, and in the application of an equivalent 
steady-flow boundary layer, these points are minor refinements to the 
theory in the experimental range considered. Nevertheless, suitable 
corrections should be applied as higher diaphragm pressures are used if 
good correlation is desired. 

The curves of figure 7^ which were comiputed on a basis of a^ = a^^ 
and Cf = O.O581R0'‘1/5^ su mma rize shock attenuation prediction by the 
linear theory for shock pressure ratios up to 10. For Cf proportional 

to this cwrve is of the type 
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where ^ is the shock-pressure attenixation at a distance l from the 
diaphragm station. According to this equation the attenuation at a 
constant nimiber of hydraulic diameters from the diaphragm would vary as 

2 

n+3 

. This relationship is dis- 
cussed more fully in a subsequent paragraph . 

Variation of pressure with time at a given distance .- Press\ire-time 
ciarves obtained experimentally by the capacitor pickup located a distance 
of 8.15 feet from the diaphragm station are compared in figure 6 with 
theoretical results computed (see fig. 8) for a distance of 8 feet and 

Cf = 0.058lRe"^/^. The error introduced by the variation in distance 

is about 1 percent of the attenuation pressxire change and may be ignored. 
The experimental shock pressvire was also obtained for each run by com- 
puting the pressure at 9*23 feet from simultaneous velocity measurements 
a r)(j then interpolating back to 8 feet along a curve parallel to the 
experimental ciirve of attenuation against distance. These velocity 
measurements fall on the theoretical measxirements within the acc\iracy of 
the plotting of figure 6. 

Both theoretical and experimental curves show a pressure increasing 
with time. This increase has also been noted when NACA miniature electri- 
cal press\n:e gages, as well as piezoelectric pressure pickups, have been 
en 5 )loyed. The initial press-ure rise across the shock determined from 
the velocity measurements, as well as from the theory, appears to fall 
below the average of the initial damped high-frequency oscillations of the 
condenser pickup pressure-time curve. However, the dynamic response of 
the Rutishauser system is not known, whereas the shock-velocity measuring 
system has a known high precision. Consequently , the use of the initial 
shock-velocity pressure point in conjunction with the Rutishauser records, 
after the high-frequency oscillations have subsided, appears a logical 
choice to give a better representation of the actual press\xre-time phe- 
nomena. Such an "experimental" curve agrees very favorably with the pres- 
sure rise predicted by the linear theory. 

Variation of other parameters with time at a given distance .- Fig- 
ures 9 to 12 show the predicted variation with time of average sonic 
speed, fluid velocity, density, and Mach number at x = 8 feet. The 
change of these parameters with time increases with shock strength, since 
the skin friction and the temperature difference between walls and fluid 
increase with shock strength for low values of shock pressure ratios. 

These aforementioned qmntities were not measiured in the experiments of 
this investigation; however, the phenomena of density increasing with 
time were found experimentally by the use of a chrono-interferometer and 
are reported in reference 8. The behavior of the average particle velocity 


the shock-tube Reynolds number 
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has not been ascertained. An increase in speed of the leading edge of the 
mixing zone which replaces the entropy discontinuity between the hot and 
cold gas regions of the theoretical flow has been reported in reference 7 
and elsewhere. However, the apparent velocity increase of the leading 
edge is also attributable to turbulent mixing since it is known that the 
width of the zone increases with time. Consequently, the velocity of the 
leading edge does not give a tnxe indication of the average particle veloc- 
ity due to the masking effect of mixing. An increase in Mach number with 
time was found experimentally in reference 5. Also, an attempt to meas- 
ure sound speed was reported in reference 5^ but the reported result of 
constant sound speed with ±2 percent variation over a 100-microsecond 
interval does not appear significant in view of the procedure employed 
and, in addition, the predicted linear-theory variation (fig. 9) is less 
than scatter of the data of reference 5« 


Correlation of Wave System and Flow Phenomena 

Wave generation and identification .- There are two ways to consider 
pressure waves: The more familiar method is that of an observer in a 

fixed position who (with waves assumed to be coming from only one direc- 
tion), upon detecting an increasing (or decreasing) pressure with time, 
recognizes the arrival of a compression (or expansion) wave. The other 
method is that of an observer, traveling with the speed of the wave, 
who (waves again assimied to be coming from only one direction) , upon 
detecting an increasing pressure with time, recognizes the generation 
of compressions. This moving observer may also determine the type of 
wave with, which he is associated by ascertaining at a given instant of 
time whether the pressure is lower (or higher) ahead of the wave, in 
which case he is traveling with a compression (or expansion) wave. 

The stationary observer will detect waves only in unsteady flow, 
whereas the moving observer will detect the growth of waves in a steady 
flow. As an example of the latter case, consider the pressure drop of 
a steady flow of air through a friction zone. It may be treated as a 
system of downstream (p) waves which, although of zero strength before 
entering the zone, become stronger expansionwise the farther into the 
zone they traverse. Simultaneously, there is a train of upstream (q) 
waves which, also of zero strength initially, grow into stronger com- 
pressions as they travel up into the friction zone. These two systems 
of moving waves combine to form a standing wave, invariant with time, 
or a pressure drop along the flow. Consequently, a fixed observer would 
recognize no moving waves in the flow. It might be noted that these 
trains of waves are the reflections of the waves which produced the 
steady flow from quiescent air. 

Cognizance should also be taken of the fact that the two observers 
would disagree on the type of supersonic Q waves in a flow. Consider 
an unsteady compression wave moving upstream against a supersonic flow 
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but actually being washed downstream relative to a fixed coordinate 
system. The traveling observer looks ahead (upstream), sees lower pres- 
STires, and so recognizes a compression wave. The fixed observer, how- 
ever, first encounters the trailing edge (high pressure) and then at a 
later instant the leading edge (low pressure) of the wave and apparently 
recognizes an expansion wave, since the pressure falls with time. Con- 
sequently, in supersonic flow the fixed observer must make additional 
velocity measurements in order to interpret the wave sign correctly. 


The imsteady waves generated in the various regions of the shock- 
tube flow can best be understood by consideration of the following equa- 
tions which result from equations (l), (2), (5)> and differentiation of 
equation (31) 
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(60) 


It is evident that the term — - is not connected with pressure or 

6t 


velocity waves but merely accoimts for the change in 


= 5 a t U due 


to traversal of the characteristic through fluid of different entropy. 

The pressure-wave generation is represented by the change in the difference 
of these terms and is a result of the skin-friction drag and heat transfer 
as expressed in the first and second terms, respectively, within the braces 
on the right-hand side of equation (60). 


If equation (60) is now examined from the traveling observer's 


viewpoint, which is more appropriate since 


sU 


represents the change 


with time as noticed by the traveling observer, the following facts 
become obvious: 

(1) Upstream (q) compression waves are generated by the effects 
of both heat addition to the fluid and skin-friction drag. 
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( 2 ) Downstream (P) compression waves are generated by heat addition. 

( 5 ) Downstream (p) expansion waves are generated by the skin-friction 

drag if M < — i — • 

7-1 

(4) Downstream (p) compression waves are generated by the skin- 

friction drag if M > — i — • 

7-1 

This last result (item (4)) was unexpected since the drag is usually 
considered as having a braking or decelerating (expansion) effect on 
the fluid ahead of it. 

The following facts mi^t be noted at this point concerning the 
heat-transfer and skin-friction effects. Expansion waves will be gen- 
erated in a for shock pressure ratios up to about 5*9 when air at 
initially uniform temperature throughout is used on both sides of the 
diaphragm. For pressure ratios above that value, compressions are gen- 
erated. Also, heat transfer is to the fl-uid in region a below shock 
pressinre ratios greater than 7-5; above this ratio the adiabatic recovery 
temperature begins to exceed the wall temperature for r = 0.9> 

For M < — - — the following conditions then apply for Too = Tg = 

7-1 

Twaii* In region a heat is being transferred to the fluid. Thus, 

the effects of heat -transfer and skin-friction drag are additive with 
regard to the generation of Q compression waves. However, in regard 
to the generation of P waves, the effects are in opposition since the 
heat transfer tends to generate compressions and the drag to generate 
expansions. In region 3 heat is being lost from the fluid to the 
walls. The heat transfer and drag are, therefore, cumialative with respect 
to the production of P expansions and in opposition with respect to 
Q waves. The stronger waves generated will, consequently, be Q compres- 
sions in region a and P expansions in region p, where the drag and 
heat-transfer effects are cumulative. 

Solution for shock pressure ratio of 2.6 . - The explanation for the 
flow behavior may now be considered with these points in mind. Bie plot 
of X against t for a shock pressiire ratio of 2.6 is shown in fig- 
\ire 15 . Adjacent to various points are indicated the values of P', Q', 

a', U, and p. The primed quantities refer to the values that the param- 
eters would have if an imaginary lamina of fluid with entropy equal to 
that of the unperturbed flow, but with velocity and pressure identical 
with the local perturbed flow, were inserted in the flow at these points. 
In other words , let 

S £ 

1 5P' _ 1 5P 1 a R 

ae 5t ag 6t 7 a^ 6t 
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and so on. Then since, by use of the imaginary lamina, the entropy term 
meikes no contribution, P' represents only the pressure and the velocity 
wave effects. 

The reason for the decrease in shock pressure with time is now evi- 
dent. Just behind the shock the overtaking P expansion waves are growing 
in strength (P' decreasing) in proportion to the time interval that it 
takes for the P characteristic to traverse the regions a and pj and, 
of coiirse, this time interval increases as regions a and p diverge 
with time. 

In figure I3, the values of P' on the p side of the entropy 
discontinuity are gradually diminishing with time from the value 
P'q = Pp = 6.^k^la„ as weak expansions are generated in region a. 

These expansions are weak because the skin-friction expansion effect is 
only slightly stronger than the opposing heat-transfer compression effect. 
The value of P' slowly decreases with time along the entropy discon- 
tinuity as the region a in which the waves are generated becomes larger. 

In region p, however, both the skin friction and heat transfer 
combine to generate strong P' expansions which increase in strength 
(P' decreasing) rapidly with distance traversed in p. Thus, the first 
P characteristic shown in figure 13 has experienced a decrement in P' 
from Pq = 6.5451 of only 0.0054 at the entropy discontinuity and of 

0.1167 when it overtakes the shock. The second P characteristic shown 
has a decrement in P' of 0.0105 at the entropy discontinuity and of 
0.2052 when it finally overtakes the shock. Note the much greater 
expansion generated in region p. 

The cause of the variation of pressure with time at a given station 
is not so obvious but may be described in a simple manner for the pax- 
ticular pressure ratio of 2.6o as follows (the explanation is not gen- 
eral since the signs of some of the waves may change with diaphragm 
pressure ratio): At a given value of xp, as time Increases the inci- 

dent downstream characteristics have traversed longer distances in a 
and shorter lengths in p. Now since the drag and heat transfer do not 
reinforce one another in a, whereas they have an acc\xmulative effect in 
P on P waves, the effect of decrease in traversal time of p more 
than offsets the effect of the increase in that of a. 

If figure 13 is again used, the following tabulation can be made for 
the values of the variation in P' along the characteristic for the four 
characteristics. Subscripts denote regions responsible for the variations. 
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Characteristics 

line 

(fig. 13 ) 

(a) 

6P*p 

Z6P' 

1 

- 0 . 005 i<- 

-0.1115 

-0.1167 

2 

-.0105 

-.0602 

-.0707 

3 

-.0179 

-.0255 

-.ol ^-54 

k 

-.0232 

-. 0047 

-.0279 


®-Value Includes reflection at entropy- 
discontinuity . 

Consequently, even though expansions are being generated in a and p, 
to a fixed observer at a given xp, the incident downstream waves appear 

as con 5 )ressions in that the value of the parameter P' is increasing 
with time simply because its reduction from the constant value Pg is 
decreasing with time. In other words, the expansion effect generated 
along the P characteristics up to a given xp is decreasing with time 
though the net expansion effect up to the shock wave is increasing 
with time. In addition, the incident Q waves in region 3 are compres- 
sion waves (for exanple, drag effect is larger than heat-transfer effect 
for this particiilar case) the strength of which increases with trav- 
ersal distance (time) from the shock to a fixed Xp. The result of 

these incident waves is thus a pressure rising with time. 

The average velocity is increasing with time at a given xp because 

> 0, downstream are stronger 

X 

than the upstream compressions. The average velocity of the entropy 
discontinui-ty decreases with time because stronger expansions (lower 
values of P ' ) are overtaking it from behind while simiiltaneously stronger 
compressions (higher values of Q') are meeting it from ahead. Both waves 
are decelerating influences. The average velocity just behind the shock 
is, of course, decreasing with time and distance. For this particular 
case, the perfect-fluid velocity is never attained in region 3 for 
t > 0. 

Althoxjgh the primed quantities of figure 13 give correct values of 
pressure and averaged velocity, the value of a' is not an indication 
of the correct average sonic speed since the primed quantities are 
always evaluated in an imaginary fl-uid lamina of constant entropy. If 
the entropy correction were applied, it would be fo-und that the average 
sonic velocity (temperature) is decreasing with time in region 3 
because the effect of the heat transferred out of the fluid is greater 
than the comibined effect of the compression waves and frictional dissi- 
pation. In other words, although the strength of the incident "compression" 


the "effective compressions," that is 
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waves increases with time at a given x, the fluid arriving at that x 
has been in motion for increasingly longer periods and more energy is 
lost by heat transfer to the walls, so that the net result is a decrease 
in temperature. 


Assessment of Various Assunrptions 

The favorable correlation between experiment and linear theory sub- 
stantiates the use of a small -perturbation approach. However, the limi- 
tations of the assumptions used in the analysis must be recognized and 
weighed carefully before application of the theory to a particular 
problem. 

One of the most precarious assumptions introduced was the averaging 
process, whereby the quasi-one-dimensional flow was treated as a pure 
one -dimens tonal flow in which averaged values of p and U were defined 
as those satisfying the continuity equation integrated over a cross sec- 
tion of the flow. It is known that this assun^jtion will Introduce errors 
in the momentum equation. These errors will be smaller for "full" veloc- 
ity profiles and, if the boundary-layer profile is assumed not to change 
with flow travel, will increase in magnitude as the ratio of cross- 
sectional average velocity to maximum velocity decreases. The errors 
will also increase as the velocity profile changes shape with flow dis- 
tance. For example, consider an incompressible steady flow of fluid 
entering a pipe with a purely one -dimensional rectangular velocity pro- 
file which later becomes a parabolic velocity profile some distance along 
the pipe. The flow has a constant average velocity, yet it requires a 
pressure drop in excess of that needed to overcome friction to account 
for the changed velocity profile. Consequently, the averaging process 
employed in deriving the basic equations may limit the application of 
the theory to boundary layers which are small relative to tube height 
and width or which have nearly constant shapes along the tube length. 

In addition, if the skin-friction integral is to be evaluated by 
the use of an equivalent-flat-plate steady flow, there is introduced 
the assumption of a shock-tube potential flow boimded by a viscous bo\xnd- 
ary layer. Even without the restrictions of an averaging process, this 
assumption requires a boundary-layer thickness 6 small in relation to 
the dimensions of the shock- tube width. This restriction was violated 
in the numerical evaluation of the theory when applied to the experi- 
mental results, since the theoretical value of 5 at the entropy dis- 
continuity was approximately 1 inch in the qi _ t>y 2-inch shock tube for 

X = 8 feet for the worst case (shock pressure ratio 1.9^). The fact 
that agreement was still good between theory and experiment even under 
these conditions indicates that this boundary assumption may be violated 
markedly without severely penalizing the accuracy of the theory. Of 
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co-urse; vhen the houndary layer develops to such a degree that it fills 
the tube and so-called "pipe flow'^ arises, then the character of the 
flov will be so altered that a new approach based on pipe flow would 
probably be more applicable than the boundary-layer and potential-flow 
approach of the equivalent-flat-plate steady flow. 

The assumption of a flow model in which velocity and thermal 
boundary-layer effects are included by averaging across the channel area 
is open to further question. This model presumes instantaneous transfer 
of heat, momentum, and so forth, from the fluid adjacent to the wall to 
that in the center of the flow. In reality, however, either molecular or 
macroscopic motion of the particles or transverse pressure waves is nec- 
essary for such a transfer so that the influence is felt at the center 
later than at the walls. Such a delay is usually of second order and 
can be ignored. 

The accuracy of the approximation for the frictional dissipation as 
equal to the product of wall shearing stress and averaged velocity cannot 
be verified unless the velocity profile is known; however, this approxi- 
mation should be good for many cases and is exact for both true one- 
dimensional flows and Poiseuille flows. 

The linearization of the differential equations prevents the appli- 
cation of the theory to large deviations from the theoretical flow. The 
degree to which this restriction may be stretched is still unknown since 
the experimental data compared with theory had variations from perfect 
fluid flow only up to fifteen percent. 

In the analysis it was assumed that the leading edge of the expan- 
sion wave traveled with a velocity -a^ and that immediately behind the 

leading edge the fluid and sonic velocities were U and a^, respec- 
tively. Since the sonic speed varies from a^ to a^ and the fluid 

velocity from 0 to U as the expansion wave is traversed, the above- 
mentioned assumption will introduce an error which will magnify as the 
expansion fan becomes larger at the higher diaphragm pressure ratios. 
Furthermore, if a more exact treatment were desired for the expansion 
region, it would be necessary to consider an equivalent accelerating 
steady flow with a favorable pressure gradient. 

In the numerical evaluation of the theory, it was assimed that the 
functional form of the skin-friction dependency on the flow length | 
was identical for the unsteady flow in the shock tube and the equivalent- 
flat-plate steady incoii5)ressible flow. This assun^tion has no rigorous 
argument. Some support may be found in the similarity of the values 
of Cf obtained from the Rayleigh, Blasius, and modified Blasius solu- 
tions. For laminar boundary layers all three solutions give functional 
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relationships, cf = (Constant) but the constants differ. No com- 

paxable solutions are known for turbulent boundary layers. Until some 
further knowledge is obtained about -unsteady boundary layers, it appears 
that this assuD 5 )tion of the equivalency of steady and vuasteady flows 
must be accepted principally on the basis of agreement between theory 
and experiment in this range. 

When the n-umerical evaluation of the theory is based on incompressi- 
ble bo-undaxy-layer theory, there is, of co-urse, a so-urce of error as dis- 
cussed earlier. However, this error can be eliminated by the use of 
closer compressible boundary-layer approximations in the evaluation of the 
Cf(|) integrals. It sho-uld be noted that the effect of the compressi- 
■bilxty correction is not as large as might first be expected. Since the 
Mach miTnh pr in region p cannot exceed I.89 for air, altho-ugh it may go 
to infinity in region a, the correction to Cf for con^jressibility will 

be much larger in a. However, Vq, becomes increasingly less than vp 
as Mq, increases, so that the influence of region a decreases. In 
addition, the compressibility effect will be less on the shock-wave 
attenuation than on the variation of the flow parameters with time at a 
given station since the iniportance of region a on 5P' changes as 
discussed in the previous section. 

The equl valent-steady-flow model does not give a_ true picture of the 
region near the entropy discontinuity even if turbulence is neglected. 
Since v is different on each side of the discontinuity, the model yields 
different boundary-layer thicknesses or an impossible discontinuous bound- 
ary layer across the discontin-uity. This effect, plus the neglected heat 
transfer across the entropy discontinuity, apparently is small when 
examined in the light of the experimental results. 


Application of Theory in Shock-Tube Design 


The compromise required between aerodynamic and mechanical design 
of shock tubes for use at high temperatures and Mach numbers becomes 
evident on scrutiny of the functional form of the equation relating the 
ratio of the attenuation shock loss or pressure rise (or fall) with time 
to the ambient pressure. This equation, for Cf proportional to 

j^-2/(n+5), is 


Poo 
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At a given value of Pp/Poo^ order to have relatively long flow dura- 

tions, the ratio Z/Oco shoirld be large; but insulating problems usually 

dictate a«j in the atmospheric range so that Z, consequently, becomes 
large. 

Then in order to minimize the timewise variations at a given x, the 
ratio z/d should be as small as possible, a requirement prescribing 
maximum available D. If D is large, mechanical strength of the shock- 
tube walls will limit the peak shock pressures so that, for lange values 
of" Pp/Poo^ i'*' will be necessary to make p„o small. A small value of p^„ 

will in turn cause to increase so that the reciprocal of the Reynolds 

number Voo/ajJ) rises and with it the attenuation and flow fluctuations. 
Consequently, the design of the shock tube will be a compromise between 
mechanical and aerodynamic design, with the z/D term (which has a more 
powerful exponent than the Reynolds number) probably being the principal 
aerodynamic consideration. 


CONCLUDING REMARKS 


Comparison of theoretical and experimental results on the flow in 
shock tubes appears to substantiate, in the range of the experiments, 
the method of a n al y sis enployed. The application of the basic concepts 
to much higher press\ire ratios and larger perturbations appears logical. 
However, in such an application, refinement of the linearized procedure 
may be required to compensate for the crudity of scane of the assumptions 
or it may even be necessary to revert to a step-by-step integration of 
the basic nonlinear differential characteristic equations and to employ 
simultaneously skin-friction coefficients based on local velocities. 

In addition, the basic characteristic equations themselves may require 
further modification as the shock strength increases to values where 
dissociation, ionization, and relaxation effects become important. Below 
the shock strengths where these last -mentioned effects occur, the trends 
predicted by the linear theory and evaluated in this paper should apply 
although their magnitude is dependent on the form of skin-friction law 
assumed. 

Although an equivalent stead.y-flow turbulent bo\indary-layer solu- 
tion gave good agreement with the experimental data, application of the 
theory to cases where the flow Reynolds number is lower than that of 
these experiments may require the use of laminar boiondary-layer solu- 
tions. Such lower flow Reynolds numbers would arise either from shocks 
weaker than those of the experiments advancing into air at atmospheric 
pressure or from shocks of the same or higher pressure ratio advancing 
into air below atmospheric pressiire. It should also be noted that, at 
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high Reynolds numbers (> 2 x 107), the evaluation of the skin-friction 

coefficient as O.O 58 I X (Reynolds number )"V5 becomes inaccurate and 
should be replaced by a more appropriate relation. 

The theoretically predicted trends in the hot gas flow for air at 
uniform temperature throughout initially are as follows : Shock pressure 

ratio decreases with distance; static pressure and average values of 
density, fluid velocity, and Mach number increase with time at a fixed 
point; and average sonic velocity decreases with time at a fixed point. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , December "J , 195^* 
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APPENDIX A 

DERIVATION OF CHARACTERISTIC EQUATIONS FOR QUASI- 
ONE-DIMENSIONAL CHANNEL FLOW WITH AREA CHANGE, 
FRICTION, AND HEAT ADDITION 


The general form of the quasi-one-dlmensional continuity equation 
is expressed as 


^pA) + ^pUA) = 0 

dt dx 


(Al) 


S loge A a loge A S loge P . ^ loge P . Bu ^ uo\ 

+ U + + U + — = O {,AZ) 

ot OX ot ox ox 


If averaged values across a section of the flow are employed for p 
and U when channel flows are considered, the area term A represents 
o nl y the physical cross-sectional area and boundary- layer displacement 
thickness is of no concern. 

The identities 


a^ = 7 RT 


(A3) 


and 


S _ d(Heat added) 
R RT 


(A4) 


combined with the equations of state and energy yield the following 


familiar relationships: 

p = (Constant) 


p = (Constant) 


2cv/R g-S/R 

(A5) 

2cp/R g-S/R 

(a6) 
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Now, if all the fluid under consideration has been in the same 
state at any previous time, equations (A5) and (a 6) may be differentiated 
with respect to x and t with no contribution from the constant term. 


S logg 

p 

2c^ 

S logg a 

. ° ^ 

Sx 


R 

Sx 

Sx 

S logg 

p 

2c^ 

S logg a 

s ^ 

° R 

St 


R 

St 

St 

S logg 

p 

ro 

o 

S logg a 

^ R 

Sx 


R 

Sx 

Sx 


(A7a) 


(A7b) 


(A7c) 


S logg p 

at 


2c a log a 
£ e 

R at 



(A7d) 


The moment-uin equation is then written by use of the skin-friction 
coefficient and hydraulic diameter as 


(a8) 

St Sx P Sx D 

This equation is not exacts since the effective values of p and U 
were chosen to satisfy the continuity equation j however, for ”full" 
velocity profiles (that is, profiles in which the velocity deviates only 
slightly from its maximum value in a large part of the cross section) , 
equation (a 8) is a good approximation. 

Equations (a 8) and (A2) can be rearranged, after substituting from 
equation (A7) and employing the convective derivative 


D( ) 8( ) , 5( ) 

St Sx 


Dt 


(A9) 
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1^5 


to give 


DU 2c^ ha. 

Dt R Sx 


^ I 2U^Cf 

+ 1 - = 0 

7 dx D 


(AlO) 


Da + a ^ - a + a ° ^ = 0 (All) 

R Dt 3x Dt Dt 


Equations (AIO) and (All) may be added and subtracted to give the 
following equations in operator form: 





a + U 


-a + a l-k + 

Dt 7 \Dt hx/R 



2Tj2c^ 

D 


(A12) 


^ \ l°Se A _ a/D S \s . 7 - 1 „ ° R , 

hxj\B. j Dt 7 \pt hxjR 7 Dt D 

(AI5) 

The assumption has been made in equations (A12) and (A13) that Cy/R 

is constant. Now, from the definition of the derivative along a charac- 
teristic of slope 6 x/ 6 t equal to U ± a as 



5 ( ) 

6 t 


dt 


+ (U ± a) 


^( ) 

dx 


D( 


Dt 


^( ) 

dx 


(Al4) 


equations (A12) and (AI 5 ) may be written in the final form: 


6 t 


l &X a * U 
5t \R / 


D log A 6 I 
e a R 

-a h + 

Dt 7 6 t 



2u2cj 

D 


(AI 5 ) 
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Iflf 


a - U 

6t 6t V R 


-a 


D logg A 
Dt 



7 6t 


+ 


7 - .J: 


7 



Dt 


2U^c^ 

~D 


(ai6) 
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APPENDIX B 

EVALUATION OF CONVECTIVE DERIVATIVE OF ENTROPY 


The convective derivative of entropy may be expressed in the form 



RTlDt 



(Bl) 


where Dj/Dt represents the heat added per unit mass by heat transfer 
and heat sources while DHp/Dt represents the heat dissipated per unit 

mass by friction. 

The heat transfer per unit mass from the walls may be approximated 
as 4h0/pD where 0 is the difference between the wall temperature and 
the recovery temperature of the flow; for example^ 


® “ ^wall ^flow 


+ r 



(B2) 


The assumption is next introduced that the relation between heat 
transfer and skin friction is given by a modified form of Reynolds ' 
analogy which is assumed to apply for turbulent as well as for laminar 
boundary layers in both steady and unsteady flow. Substitution of the 
perfect gas law into the modified form of Reynolds ' analogy resiilts in 
the following equation for the wall heat-transfer coefficient: 




The heat-transfer term may then be written as 


(B3) 


Dt 


27 ^ 

7 - 1 D 



(Bk) 
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The frictional dissipation per unit mass is assumed to be the product 
of the average velocity mioltiplied by the wall-shearing force and divided 
by the mass in the flow: 


^2f _ Ut( P erimeter) Ax _ 2U^ 
Dt pA Ax D 


Cf 


(B5) 


Although equation (B5) is not exact in general, it is a good approxi- 
mation for most cases. It might be noted that equation (B5) is exact for 
a true one -dimensional-flow model where the velocity is considered con- 
stant across a cross section of the flow with a discontinuity in velocity 
at the wall. It is also exact for incompressible Poiseuille flow since 


the dissipation function, 
tion (B5) • 




2 

I dy, reduces to the value of equa- 


Therefore, under the conditions presumed to apply above, equa- 
tions (B1), (b 4), and (B5) may be combined to evaluate the convective 
derivative 




^ i 

7 - 1 T y D 


(b6) 
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APPENDIX C 

LINEARIZED SOLUTION FOR INTERSECTION OF CHARACTERISTICS 
AND ENTROPY DISCONTINUITY 


The linearized solution for the intersection of the characteristic 
lines with an entropy discontinuity is obtained in the following manner 
(see fig. 1 for notation); 


6P = Pf - Pc = Qf - Qc = - 1^ (Cl) 

Pc + Qf .= -j— + ij (C2) 


In these equations the values of 7 are assumed to be Identical 
everywhere . 

Dividing equation (Cl) by (C2) and applying the condition that the 
pressures are always equal to one another at any instant on both sides 
of the discontinuity yields 


6P 

Pc + Qf 



Sf “Sc 



(C3) 


Now, if the subscript o denotes conditions at time t = 0 immedi- 
ately after the diaphragm burst with 


^ = Sf - Sc (C4) 

or 

(C5) 


ZiS = 2So(l + A) where A« 1 


1+8 
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(c6) 


substituting equation (C5) in. equation (C3) and expanding gives 



5P 

5Po 


(Pq - + (Q-f ~ ^Po) 

^cco + 


1 + A^ + 



(C9) 
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Pq - Pet Qjp - Qpq 

Since and are of order 

Pqo + '^Po ^oo + ^Po 

or linearized solution becomes 


A, the first-order 


6P _ 6Q 

8P^~ 5Qo 


(Pe - P«o) ^ («f - «So) 

1 + 1 1 L + 

^oo + ^Po 


7 - 1 

7 






(CIO) 
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APPENDIX D 

RELATION BETWEEN SKIN-FRICTION COEFFICIENT 
AND VELOCITY PROFILE 


The relation between skin-friction coefficient and velocity profile 
for an incompressible turbulent boundary layer on a flat plate may be 
found in the following manner. Let the velocity profile in the boundary 

^ / \l/n 

layer be of the form -2L = ( Y. ) where u^ is the local streamwise 

U \6/ ^ 

velocity at the distance y from the plate and 6 is the boundary- layer 
thickness. Then, equating the wall shearing stress t to the rate of 
decrease, in the streamwise | direction, of the momentum in the boundary 
layer gives 


T = — / pUy(U - Uy)dy 

d| Jq 

Evaluation of this equation for constant p yields 

T _ n ^ ^ (2 - n2)s dn 

pU^ (n + l)(n + 2) d| (n + i)2(n + 2 )^ 


(Dl) 


(D2) 


The following relations are then assumed between the shearing-stress 
velocity U*’ and the boundary-layer thickness 6 (see ref. 10): 



(D5) 




i, 


The parameter 
but its form for 


B may be assimied a constant for 
n = n(|) is not specified. 


(d4) 


n invariant vith 
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Subs'tl'tu'tion of equations (D3) and (D^) into equation (D2) results 
in the following differential equation; 


gn+1 d.5 


1 Uv 

\ B 


l/nl 


2n 

,n+l 


(n + l)(n + 2) 


o n+3 

(2 - n2) gn+1 ^ 


n 


n(n + 1) (n + 2) 


(d5) 

For the case of n and B independent of equation (D5) may be 
Integrated, by assming n -1. This integration with application of 
the bo\mdary condition of 6=0 at | = 0 produces 


n+1 


6 - — 

(n + 2)(n + 3) 

UB^ 

n 

Substitution of equation (d 6) into 


n+3 


n+1 

n+3 


(d6) 


expression for the skin-friction coefficient: 


puy2 


n 


(n + 2)(n + 3)B^ 


n+3 


2 

n+5 


(D7) 
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Figure 1.- Distance-time plot of shock-tube flow. 
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Figure 2.- Schematic diagram of test equipment. (All power voltage regu- 
lated at 115 volts, alternating current.) 


MCA TO 3575 



L-85316*1 


Figure 5.- Photograph of measuring station in shock tubes. 
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Amplifier (2D2I) 
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reset switch 



-O + 250 V 
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^ Anode 
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R| 22K 
RZ 22 K 
R3 22 K 
R4 100 K 
R5 100 K 1/2 w 


Re 22K 2w 
Ry 470 K l/2w 
R0 15 K 2 w 
R9 1.2 megohm 
R|0 3.3 K 2 w 


1/2 w 


O Output to 
counter 


C| 300 /i./xfd 
Cz 300 
C3 0.001 /xfd 
C4 16 /xfd 450v 
C5 0.001 /i.fd 600 V 


Figure 4.- Circuit diagram of phototube and amplifiers. 
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(a) Pg = 45 Ib/sq in. gage. Experimental data corrected to 

Pg/Pco “ 4.o6l. 

Figure Theoretical and experimental shock pressure ratio plotted 

against distance from diaphragm. 
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(b) Pg = TO IP/sq in. gage. Experimental data corrected to 

Pg/Poo = 


Figure 5*- Continued. 
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(c) Pg = 95 Ib/sq in. gage. Experimental data corrected to 

Pe/Pco = 7^55.- 


Figxire 5*- Continued 
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(d) Pg = 250 Ib/sq in. gage. Experimental data corrected to 

Pe/Pco = 17.915. 


Figure 5 .- Concluded 
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(b) Pg /p„= 5.764, 




Time 


(d) Pg /P^= 17.915. 

Figure 6.- Con 5 )arison of experimental and theoretical pressure-time curves. 
Xgxp = 8.13 feet; = 8.0 feet. Circles represent theoretical 

values. The sine-wave time base has a frequency of 1^000 cps. 
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(a) Shock pressure -loss coefficient. 

Figure 7 .- Variation of shock pressure-loss and shock pressure-attenuation 
coefficients with perfect-fluid shock pressure ratio, = aooj n = 7J 

Cf = 0.058lRe"^/5. 
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(b) Shock press\are-attenmtion coefficient. 


Figure 7-- Concluded 
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Figure 8.- Theoretical variation of static pressure ratio with time. 
X = 8 feet; a^ = aoo = 1A17 feet per second; Cf = 0.058lRe“^/^. 
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Figure 9 *- Theoretical variation of average sonic -velocity ratio ay/S’e 

With time, x = 8 feet; = a^o = ^*eet per second; Cf = 0.058lEe“'^/^ . 
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Figure 10.- Theoretical variation of average particle- velocity ratio Uv/^e 
with time, x = 8 feet; = a^ = 1,11? feet per second; Cf = 0.058lEg"^/5. 


I 


MCA TN 5375 


Average density ratio 



Figure !!•- Theoretical variation of average density ratio with time. 
X = 8 feet; a^ = a^^^ = 1,117 feet per second; Cf = 0 . 058 lR^“’^/ 5 . 
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Figure 12.- Theoretical variation of average Mach number with time. 
X = 8 feet; = a^o = 1,117 feet per second; Cf = 0.058lFe'"^/5. 
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Figure IJ.- Characteristic plot for shock pressure ratio of 2.6o showing 

generation of waves in region p. 
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